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ABSTRACT 



o 



Context. Planets are supposed to form in circumstellar disks. The additional gravitational potential of a planet perturbs the disk and 
leads to characteristic structures, i.e. spiral waves and gaps, in the disk's density profile. 
Aims. We perform a large-scale parameter study of the observability of these planet-induced structures in circumstellar disks in the 
D ' (sub)mm wavelength range for the Atacama Large (Sub)Millimeter Array (ALMA). 

, Methods. On the basis of hydrodynamical and magneto-hydrodynamical simulations of star-disk-planet models, we calculated the 
disk temperature structure and (sub)mm images of these systems. These were used to derive simulated ALMA images. Because 
appropriate objects are frequent in the Taurus- Auriga region, we focused on a distance of 140 pc and a declination of « 20°. The 
t-H ' explored range of star-disk-planet configurations consists of six hydrodynamical simulations (including magnetic fields and different 
planet masses), nine disk sizes with outer radii ranging from 9 AU to 225 AU, 15 total disk masses in the range between 2.67 • 10~ 7 M Q 
r ^" 1 . and 4.10 • 10~ 2 M Q , six different central stars, and two different grain size distributions, resulting in 10 000 disk models. 

Results. On almost all scales and in particular down to a scale of a few AU, ALMA is able to trace disk structures induced by 
planet-disk interaction or by the influence of magnetic fields on the wavelength range between 0.4 and 2.0 mm. In most cases, the 
optimum angular resolution is limited by the sensitivity of ALMA. However, within the range of typical masses of protoplanetary 
disks (0.1 M©- 0.001 M Q ) the disk mass has a minor impact on the observability. It is possible to resolve disks down to 2.67 • 10~ 6 M Q 
and trace gaps induced by a planet with jf- = 0.001 in disks with 2.67 • 10" M Q with a signal-to-noise ratio greater than three. 
The central star has a major impact on the observability of gaps, as well as the considered maximum grainsize of the dust in the 
disk. In general, it is more likely to trace planet-induced gaps in our magnetohydrodynamical disk models, because gaps are wider 
in the presence of magnetic fields. We also find that zonal flows resulting from magneto-rotational instability (MRI) create gap-like 
structures in the disk's re-emission radiation, which are observable with ALMA. 

Conclusions. Through the unprecedented resolution and sensitivity of ALMA in the (sub)mm wavelength range, the expected detailed 
^ ■ observations of planet-disk interaction and global disk structures will deepen our understanding of the planet formation and disk 
evolution process. 
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. |1. Introduction side of the orbiting planet (e.g. iPapaloizou & Li n 1984; Klahr 

Cn( . . & Kley 2006). The inner spiral propagates inwards, while the 

Although an early configuration of the Atacama Large (Sub)- outer spiral propagates outwa rds through the disk. In a frame of 
. Millimeter Array (ALMA) has become only recently available, refer ence comoving with the planet, these structures are static 
ryH us capabilities and its potential for groundbreaking dtskovenes and moye with k These spirals haye low amplitude in den . 
> have been demonstrated to be enormous (e.g.| Boley et al. ||2012hJ sity enhancement , typically < 10% of the background disk col- 
L> Already in the ALMA cycle confi g urations its an gular resolu- umn density The other pi an et- induced structure is a de nsity gap 
K tion and sensitivity is several times higher than those of the Sub- around the orbit of the planet (Goldreich & Tremaine 1980 Pa- 
millimeter Array. After completion of the entire array, ALMA paloizou & Lin 1984) . For massive planetS; the tidal torques on 
>» .will allow one to observe the density structure of young circum- the surrounding gas overcome the pressure gradient and the vis- 
stellar disks in unprecedented detail. cous spread j ng; and th e gas around the orb i tal position is pushed 
Planets are supposed to form in these disks (e.g. Mordasini away from corotation (Br vden et al. 1999 C rida etal.l 20061 de 
et al. 2010) The two main scenari os are the core accretion model Val-Borro et al. 20061). The density in the gap can be a very 
or the nucleated instability (e.g. |Mordasini et al| 20101) and the small f raction f the initial density. Even in more evolved sys- 
planet formation by a gravitational disk instability (e.g. Durisen tems like deb ris disks, giant p lanets induce characteristic disk 
et al. 2007) Planets orbiting in circumstel lar disks interact gray- structures (see lErtele t al. 2012 for a detailed study). 
itationally with the surroundi ng disk gas dGoldreich & Tremaind 

1980l; lPapaloizou & Lin 1984). Planets with planet-to-star mass Grid-based hydrodyna mical or smoothed particle hydrody- 
ratios over 0.0001 cr eate significant density perturbations in the namical (SPH) simulations dD'Angelo e t al. 2002i: lFouchet et all 
disks dTanaka et alJ e.g. 20021; IWardl 19971 IPapaloizou & Linl 20101), comb ined with follow-up radiative tran sfer simulations 
19841). One of these structures is a spiral density wave on each (IWolf e t al. 2002 ; Jang-Condell & Sasselov 200 31: Wolf & 

D'Angelo 2005 Jang-Condell & Boss 2007 ; Jang-Condell 2009 ; 



ruge@astrophysik.uni-kiel.de iGonzalez et al.1 20101; iGonzalez et al.1 2012 ), have shown that 
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large-scale disk structures are visible through the continuum re- 
emission and scattered light of the disk. The case of the ob- 
servability of simulated line emission has b een i nvestigated by 
ISemenov eta l. (2008), Re galv et all (20101), and lCleeves et al] 
(20111) For the infrared wavelength range, the observability 
of planet-induced, la r ge-scale disk struc t ures w as explored by 
[Vamiere et al. (2006), Wol f & D'Angelolf 20051) and Jang-Con- 
dell (2009) However, due to the much lower optical depth in 
the (sub)mm regime, observations in this wavelength range are 
better suited to directly tracing the density structure of the disk 
interior. But, for high angular resolution observations in this 
wavelength range, large interferometer arrays, like ALMA, are 
required. 

The intrinsic implementation of the stellar radiation field in 
the hydrodynamical simulations affects the gap form and en- 
large s the vertical extension especi ally of the outer edge of the 
gap dJang-Condell & Turned l2012l)l Furthermore, it has been 
found that in SPH-models gaps are enl arged in comparison t o 
grid-based hydrodynamical simulations (iGonzalez et al.l 2010) 
Also, gaps can be enla rged through the detachment of la rger dust 
grains from the gas (iPaardekooper & Mellemal 20041) There- 
fore, we consider non radiative hydrodynamical simulations as 
a simpler disk approximation. The additional presence of mag- 
netic fields within the disk influenc es the planet formation and 
planet-disk interaction process (e.g. Papaloizou & Nelson 2003 ; 
Nelson & Papaloizou 2003; Papaloizou et al. 2004; Nelson & 
Papaloizou 20041; iJohansen et al.l 2009; Pinil la et al] 20121). In 
circumstellar accretion disks, it is well known that molecular 
viscosity is not sufficient for the outward transport of angular 
momentum and for the transport of mass i nwards, as is required 
for mass to accrete onto the central star dPringlel 1981 Balbus 
& Hawley 1998). It was shown that hydrodynamical turbulence 
from a non linear shear-instability cannot transport angular mo- 
mentum with t he requi red efficiency to support accretion at the 
observed rates (|Ji et a l. 2006) Thus, many mechanisms have 
been proposed to explain the source of the dissipation that drives 
accretion. Large-scale magnetic fields with open field line con- 
figurations can play a role in th e vertical transport of a ngular 
momentum away from the disk (iKonigl & Pudrit z 2000) , The 
baroclinic instability has also been proposed as a possible mech- 
anism for an gular m omentum transp ort via vortices ( Klahr & 
Bodenheimer 2003) One of the most promising mechanisms for 
driving accretion, and certainly the easiest to incorporate in nu- 
merical simulations, is the magneto-rotational instability (MRI). 
Circumstellar Keplerian disks that are permeated by a weak 
magnetic field are linearly unstable (iBalbus & Hawleyl 1991) 
This instability results in the development of turbulence, which 
can transport mass and angular momentum efficiently (Balbus & 
Hawley 1998) , A minimum degree of ionization leading to good 
coupling between the gas and the magnetic field in the disk is 
required for the MRI to function. In the cold and dusty midplane 
regions the degree of ionization will be low, possibly produc- 
ing a dead zone, while in the outer parts of the disk, cosmic 
rays can penetrate, allowing for the necessary ionization degre e 
dTurner et al.l 20101; iDzvurkevich et al.l 20101 iFlock et al.l 201 II), 
yet here the low densities will make ambipolar diffusion a power- 
ful way to suppress the magneto-rotational instabilities. For the 
need to incorporate turbulence, we assume that our disks follow 
the laws of ideal magnetohydrodynamics. More realistic models 
of non ideal magnetohydrodynamical disk simulations are still 
un der development. 

IWolf et a l. ( 2002) show that ALMA is indeed able to trace 
pre planetary and planet-induced large-scale structures in cir- 
cumstellar disks and even allows one to observe the circum- 



planeta rv accretion region (IWolf & D'Ang elo 2005) Cossins 
et al. ( 2010) find that spiral structures in compact self-gravitating 
disks, which are the same size as the planet-induced large-scale 
structures, are visible in ALMA observations over a wider range 
of wavelength. While most of the earlier studies of the observ- 
ability of planet-induced large-scale structures were mainly fo- 
cused on exemplary case studies or on the impact of high plan- 
etary masses (5 MTi T . lGonzalez et al . 2012), we now investigate 
a large parameter space of stellar, planetary, and disk parame- 
ters, to exactly evaluate under which conditions ALMA will al- 
low one to trace characteristic density structures resulting from 
planet-disk interaction or magneto-rotational instability (MRI) 
best. Our model setup contains six 3D-hydrodynamical simu- 
lations (including magnetic fields and in particular lower planet 
masses), where gas and dust are homogeneously mixed. It also 
includes nine disk sizes with outer radii from 9 AU to 225 AU, 
15 total disk masses in the range between 2.67 • 10" 7 M o and 
4.10 • 10" 2 M , four main sequence and two pre-main- sequence 
stars and two different grain size distributions. Besides this, 
14 ALMA configurations, seven observing wavelengths in the 
range from 330 jim to 3.3 mm and three exposure times are con- 
sidered. 

This work is structured as follows. In § [2] we outline the 
hydrodynamical and radiative transfer simulation techniques and 
describe the parameter space of our investigations in detail in § 
[3] The data analysis is outlined in § |4] Our results are shown in 
§ [5] Finally, we present our conclusions in § [6] 

2. Simulation techniques 

Motivated bv llWolf et a l. (2002) and I Wolf & D' An geld (20051), 
who investigated the observability of planet-induced structures 
in a selected sample of circumstellar disks, we now explore 
which properties of the central star, the disk, and the planet con- 
strain clear detections. In the following, we briefly outline the 
applied hydrodynamical and the radiative transfer techniques. 
Thereafter, the observing simulations are diskussed briefly. 

2.1. Hydrodynamical simulation 

To simulate the interaction between a planet and a circumstel- 
lar disk, three-dimensional hydro- and magneto-hydrodynamical 
simulations were car ried out using the fin ite volume fluid dy- 
namics code PLUTO dMignone et al. 2007) , This computes Rie- 
mann fluxes using the HLLC and HLLD solvers for the hydro- 
dynamic and magneto-hydrodynamic cases, respectively. Time 
integration is done using a second-order Runge Kutta scheme, 
and spatial interpolation is done using a second-order total varia- 
tion diminishing (TVD) scheme. For the simulations that include 
magnetic fields, the constrained trans port method for preserv- 
ing a divergence-free magnetic field dGardiner & Ston e 2005) 
is used. Adiabatic index y = 1 and isothermal approach is as- 
sumed, e.g., T = Tq (rsinft)~ 2b see definition of speed of sound 
below. 

The geometry of the computational grid is spherical; i.e., the 
coordinates are given as (r, ft, (p). The domain in code units are 
given by re [1, 10], ft e [tt/2 - 0.3,tt/2 + 0.3] and <p e [0,2tt]. 
The grid resolution is (N^N^Ny) = (256, 128,256), and it is 
centered on the center of mass of the planet-star system. 

The gas disk is initially in sub-Keplerian rotation around the 
central star. The azimuthal velocity is given by 

v<p = yjv*-c*(a- 2b), (1) 
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Table 1. Table of all simulations. 



Name 


Simulation 


(N t ,N»,N v ) 


7] = M p /M* 




r domain 


& domain (d&) 


tp domain 


h/r 


SF 






[AU] 


[AU] 


[rad] 


[rad] 






SI 


HD 


(256,128,256) 


no planet 




[1,10] 


0.3 


In 


0.07 


unperturbed 


SP4 


HD 


(256,128,256) 


10" 4 


5.0 


[1,10] 


0.3 


In 


0.07 


gap 


SP1 


HD 


(256,128,256) 


10" 3 


5.0 


[1,10] 


0.3 


In 


0.07 


gap 


ST1 


MHD 


(256,128,256) 


no planet 




[1,10] 


0.3 


In 


0.07 


zonal flow 


SPT4 


MHD 


(256,128,256) 


10" 4 


5.0 


[1,10] 


0.3 


In 


0.07 


P-gap 


SPT6 


MHD 


(256,128,256) 


io- 3 


5.0 


[1,10] 


0.3 


In 


0.07 


gap 



Notes. In the last column, SF refers to "Special Feature" and the terms have the following meaning, p.gap means partial gap, gap means full gap 
and zonal flow means that there are magnetic field induced density bumps in the disk. The hydrodynamical simulations are characterized as HD 
and the magnetohydrodynamical ones as MHD. The name gives the link to the density maps shown in Fig. [T] 



where is the Keplerian velocity and a = 3/2 and b = 0.5 are 
the exponents of the radial power law distribution of the density 
p oc r~ a and sound speed c s = cq (rsin&)~ b . The initial density 
distribution is given by 

p(r, 0) = (r sin #)" 3/2 exp Sm ^ 2 ~ 1 . (2) 

v c o / 

The disk is described by a locally isothermal equation of state 
P = c\p. The ratio of the pressure scale height h to the ra- 
dial coordinate of the disk is taken to be a constant such that 
h = H/(r sin#) = 0.07. Therefore, the disk models are not 
flared. However, since flaring becomes important for disks pre- 
dominantly heated by the central star, flaring will be considered 
in future studies. 

The gravitational potential of the planet is given by a soft- 
ened point-mass potential 

GM 

°P (r) = ~7i 12~ 2M/2 (3) 

(Ir-r^l 2 + 6 2 ) 1 / 2 

where e is the softening parameter, needed to avoid very high 
values of the potential near the planet position that could induce 
very high velocities and numerical problems. For all the simula- 
tions e is set to be a fraction of the Hill radius e = I • r p (M /? /3) 1/3 
with / = 0.3. 

The vertical extent of the simulations is # e [7r/2-0.3,7r/2 + 
0.3], since for the purpose of the hydrodynamical evolution, 
there is no need to simulate the entire spherical domain # e 
[0,7r], because this would only complicate the simulations and 
significantly increase the computational time. The simulated do- 
main covers four pressure scale heights above and below the 
midplane of the disk. For the purpose of the radiative transfer 
analysis, however, it is necessary to have the full # domain, so 
as to follow the radiation at all possible angles. For this reason, 
the domain is increase d a-p osteriori, providing the full n range in 
the # direction (see § 12.21) . In the extended domain, the density 
in the disk is set to zero, which is a good approximation, since 
the density at four pressure scale heights has decreased several 
orders of magnitude over the midplane density. 

Additionally, since in the magnetohydrodynamical simula- 
tions, the zones r < 2 and 9 < r < 10 are boundary buffer 
zones where resistivity is applied, these parts of the domain are 
generally not included in the radiation transport analysis. This 
means that the effective grid for the radiative transfer is r e [2, 9] 
(in PLUTO code units). The number of cells in PLUTO is 
(N^NfoNy) = (200, 128,256). 

Different configurations were explored. The planet mass M p 
to stellar mass M* ratio is denoted as 

M p 

77=—. (4) 



A list of the simulations performed and of the parameters used 
in each of them is given in Table [T] Our setup contains six simu- 
lations in all. First we performed a simulation of an unperturbed 
disk to compare the other simulations with. Then we focused on 
planets with 77 = 0.001 to investigate the observability of Jovian- 
mass objects in circumstellar disks. Furthermore, we selected 
planets with n = 0.0001 to estimate the possibility of tracing 
lower mass objects with ALMA; i.e., the mass of a planet de- 
pends on the mass of the star. Additionally, all simulations were 
carried out for the hydro- and the magnetohydrodynamical case 
to determine the influence of magnetic fields on the observabil- 
ity of large-scale structures. Local and global magnetohydrody- 
namical simulations of accretion disks without the presence of a 
planet have shown the formation of zonal flow features, which 
are axis symmetric density enhancements in the disk confirmed 
in modifications of the rotation law (see Fig. [T]top, left). These 
features are effects of the non locality of the Maxwell stresses 
associated with MRI. These zonal flows are assumed to be of 
major importance for the accu mulation of dust: 1. ) retain parti- 
cles for a long time in the disk (iPinilla et al . 2012), 2.) to create 
large enough gas-to-dust ratios to trigger planetesimals forma- 
tion dJohansen et al.l 20091) , 

2.2. Radiative transfer 

Continuum radiative transfer: The follow-up radiative trans- 
fer simulations were performed with th e Monte Carlo 3D con- 
tinuum radiative transfer code MC3D (Wolf et al. 1999; Wolf 
2003), which self-consistently calculates the dust temperature 
of an arbitrary density distribution. Based on this, the result- 
ing spectral density distributions (SED), re-emission, and scat- 
tered light maps were simulated. Besides analytically defined 
dust density profiles, MC3D allows reprocessing a numerically 
generated density structure, such as those resulting from the hy- 
drodynamical simulations. However, the outcome of the simula- 
tions with PLUTO has to be converted for the use in MC3D. 

Conversion and input of hydrodynamical data in MC3D: By 

merging eight neighboring cells from the PLUTO grid into one 
we significantly decrease the required computing time, which is 
essential given the large parameter space (see § |3J. The cor- 
responding diskrepancies between a high resolution re-emission 
image and a merged one are lower than 0.1%. Additionally, the 
PLUTO code computes its resul ts in a cut-out in the ^-direction 
around the disk midplane (see § 12.11) . This is reproduced by as- 
suming a density equal to zero in the p art o f the model space 
that was not simulated with PLUTO (§ 12.1b . Furthermore we 
use as many cells in ^-direction as needed to get the same res- 
olution in the cut-out area. The grid structure of MC3D is 
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-0.90 
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-0.63 



-0.85 



-1.07 



-1.29 
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0.30 



-0.20 



-0.71 



-1.21 



-1.71 



-2.21 



-2.71 



Fig. 1. Midplane density of all hydro- and magnetohydrodynamical disk models. The units of the axis are AU, and the midplane density is given 
in logarithmic scale of the code density unit. See Table [T]for a detailed model description. 

(N T , N&, N<p) = (100, 341, 128) resulting in a total cell number Dust properties: It is assumed that the gas and the dust within 
of 4.36 • 10 6 . The images simulated in the radiative transfer are the circumstellar disk are perfectly mixed with a gas-to-dust 
fixed to a resolution of 201 px x 201 px. This ensures that the mass ratio of 100 : 1. The dust grains are spherical and consist of 
synthesized beam of ALMA in the simulated observations al- 62.5 % si licate and 37.5 % graphite (optical data bv Weingartner 
ways contains several tens of pixels (see § [23}. & Draine 2001 and concept by|MathisetaL| 1911). Ice were not 
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Fig. 2. Absorption efficiency for the large and small dust grains in the 
mid-infrared to mm-regime. The solid line represents large dust grains 
and the dashed line small dust grains. For furth er informati on see § 12.21 
All plots were created with python matplotlib ( Hunter 2007) 



Spatial resolution of ALMA 
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Fig. 3. Resolution of ALMA at the seven wavelengths considered. 
Different line styles represent different wavelengths. These values are 
derived from the synthesized beams simulated with CASA. 
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Impact of phase noise 
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25 



Fig. 4. Total-to-thermal noise ratio. For the phase noise calculation 
we use a PWV= 1.0 mm with a deviation of APWV = 0.15 mm, with 
an exposure time of two hours. Different line styles represent different 
wavelengths. Thanks to the sequential simulation process of thermal 
and phase noise, ratios below 1.00 are possible, because the phase and 
thermal noise can erase each other. 



[2] The grain size distribution can be approximated by one effec- 
tive dust grain radius a e ff resulting from the volume averaging: 



4-q _ 4-q 
a max a m [ n 



- a 



V 



(5) 



Based on this, the mass density from the hydro- and magnetohy- 
drodynamical simulations is transformed into a number density 
as needed by the MC3D. The planetary heati ng of the circum- 
planet ary environment, considered already by Wol f & D' An gelo 
(2005), and heating processes by viscosity and accretion onto 
the star or the planet are not taken into account in this study. 
The only heating source is the photosphere of the central star in 
a blackbody approximation, therefore the disks are heated pas- 
sively. 



considered in the current the study, but it will be regarded on our 
future investigations because its impact in particular in the outer 
parts of circumstellar disks is assumed not to be negligible. The 
grain size distributi on follows a pow er law, n(a) oc a~ q with an 
exponent q = 2.5 (iDohnanyil 19691), We selected two ranges 
for the dust grain radius a. The range (a e [0.005 /mi, 100 fim]) 
mimics the larger grain size frequently found in th e denser re- 
gions around the midplane of a circumstellar disk (iWolf et al.l 
2003; Madlener et al. 2012). Since giant planets induce struc- 
tures, in particular in this region of a disk, it is necessary to 
consider this dust species in our study. The other dust grain ra- 
dius rang e (a e [0.005 /im, 0.25 /im]) represents the interstellar 
medium (iSavage & Mathisl 1979), which is important for this 
study becaus e some young stellar objects (e.g. CB26. Sauter 
et al. 2009) show a lack of large dust particles. 

The absorption properties of the dust depend on the grain 
size distribution. For comparison, the absorption efficiency g a bs 
calculated using Mie theory for both grain sizes is shown in Fig. 



2.3. (Sub)mm-observations: Simulation setup 

ALMA: After its completion ALMA will consist of 66 
( sub )mm- antennas (54 x 12-m- and 12 x 7-m-antennas), which 
can be arranged in 28 different array configurations. ALMA will 
be able to observe in ten bands between 0.3 mm and 9.3 mm 
dBrown et al.1 20041) 

CASA: We simulated observations of a circumstellar disk for 
every second A LMA-antenna configuration, by using the CASA 
3.2 simulator dPetrv for the CASA Development Team! 20121) 
The number of considered bands/wavelengths amounts to seven 
(shown in Tab. O. The angular resolution of ALMA is shown 
in Fig. [3] as a function of wavelength and array configuration. 
See appendix [A| for further information on the interferometer. 
Furthermore, we consider three different exposure times (^h, 2h, 
8h). 

The simulated observations include the significant influence 
of thermal noise, which is caused by precipitable water va- 
por (PWV) in the Earth's atmosphere. For the wavelength- 
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Table 2. Selected wavelengths for the simulated ALMA observations. 



Central wavelength X 
[jim] 


Central frequency v 
[GHz] 


ALMA band 


330 


910 


10 


430 


700 


9 


700 


430 


8 


950 


320 


7 


1300 


270 


6 


2000 


150 


4 


3300 


90 


3 



Notes. Bandwidth: 8 GHz 



depende nce of the PWV we follow the CASA cookbook (Ott 
& Kern 201 ID and consider the recommended values. 

From simulated observations of an empty patch of the sky, 
noise maps are derived to obtain the impact of thermal noise on 
our setup. In addition, the CASA simulator is also able to sim- 
ulate phase noise. Because of the high computing time require- 
ment of phase noise simulations, we first explored the impact 
of phase noise on the image quality by adding it to the already 
calculated thermal noise maps (we call the outcome total noise). 
The total-to-thermal noise ratio is plotted in Fig. [4] It only shows 
significant deviations (> 10% for every ALMA configuration) 
for the shortest wavelength. 

In another test, 10% of our parameter space was simulated 
again, considering thermal and phase noise. In comparison to 
the case with thermal noise only, just 0.6% of the simulated ob- 
servations are evaluated differently (for assessment criteria see § 
[4]). For this reason we decide to consider only thermal noise for 
our study. 

Noise maps and all simulated observations were calculated 
for the Taurus- Auriga star forming region considering the posi- 
tion of TTAU (a = 04h22m59s, 6 =+19°32'06", J2000), which 
corresponds to a maximum elongation of « 50° in December. 
On the smallest scales of the spatial resolution of ALMA it can 
be necessary to figure the observation out with mosaic imaging if 
the maximum angular scale of ALMA is smaller than the object. 



3. Model setup 

Disk size : The self- similar hydro- and magnetohydrodynamical 
simulations allow us to arbitrarily scale the size and the total 
mass of the disk. The disk size is adapted by scaling the inner 
and outer disk radius linearly with a factor k • AU in the range 
ofke {1,4,7,10,13,16,19,22,25}. In the case of k = 1, the 
hydrodynamical simulations are used unsealed. The maximum 
value k = 25 was chosen to simulate a pl anet with a semi- major 
axis (125 AU) similar to Fomalhaut b (iKalas et a l. 2008), For 
illustration the inner and outer disk radius and the semi-major 
axis of the planet are shown as a function of k in Fig. [5j In a 
more general view, our disk models are anulli or smaller cut-outs 
from large circumstellar disks. However, since we investigate 
the ability to trace structures of a defined size in our disk models, 
it is not necessary to take more disk areas into account. In the 
following, we call these models disks. 

Total disk masses: The typical estimates of the total mass of 
circumstellar disks is in the range of 0.01 to 0.1 ( An- 
drews et al. 20101) During the disk evolution, the disk mass 
decreases because of, i.e., photoevaporation and accretion (Wolf 




- - Inner radius, r in (fc=l) = 2 AU 

— Semi major axis, r pl (k = l) = 5 AU 
Outer radius, r ou (k = l) = 9 AU 



100 150 
Radius [AU] 

Fig. 5. Correlation between scaling factor k and the inner (dashed 
line) and outer disk (dotted line) radius (r in , r out ) and the semi-major 
axis (kolid line) of the planet (r p i). 



et al. 20 121) Half of all ci rcumstellar disks have a lifetime of 
< 3 Myr (lHaisch et al.l 2001) , To consider this factor in our cur- 
rent study, a wide range of masses was taken into account. At the 
lowest disk masses in our setup, which mimic the mass typical 
of a debris disk, our hydro- and magnetohydrodynamical disk 
models may no longer correcty approximate the disk structure, 
because additional processes, e.g. radiation pressure and photoe- 
vaporation, become important. Nevertheless, observations were 
simulated to predict the minimum mass that allows one a de- 
tection of a circumstellar disk. The following disk masses are 
considered (only integers in the exponent): 



M = 2.67 • 10" 



(6) 



The value 2.67 results from the volume integration of the density 
structure of the disk models from the hydro- and magnetohydro- 
dynamical simulations. As an alternative to this approach, we 
prepared a setup for a total disk m ass d erive d from the disk of 
the Butterfly Star ( Wolf et al. 2003 ) in § IP 



Distance, stars, and disk inclination: Throughout the entire 
study objects are located at a distance of 140 pc (typical for TT- 
objects). We assume disks that are seen almost face-on (inclina- 
tion angle i = 5°). The number of different central stars, charac- 
terized by luminosity and effective temperature, amounts to six. 
Table [3] summarizes the important stellar parameters for the four 
main- sequence stars and two typical pre-main- sequence stars. 

Table 3. List of the selected stars. 



Spectral type 


L 


T e ff 


M 




[L ] 


[K] 


[M Q ] 


K(MS) 


0.35 


4500 


0.7 


G(MS) 


1 


6000 


1.0 


F(MS) 


7.5 


6900 


1.8 


A (MS) 


20 


8500 


2.4 


T-Tauri (PMS) 


0.95 


4000 


0.5 


Herbig Ae (PMS) 


43 


9500 


2.5 



Notes. The denotation and the effective temperature of the main- 
sequence stars as spectral types refers to the Harvard Classification of 
Stars. We use the mass-luminosity relationship to estimate the stellar 
mass. (MS) = main- sequence stars; (PMS) = pre-main- sequence stars. 
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Fig. 6. Illustration of the radial brightness distribution (solid line) with 
auxiliary line (dotted line) used for the automatical gap detection. The 
corresponding simulated observation is shown in the upper right. The 
disk model has a scaling factor of k = 22 (outer disk radius 198 AU), 
a total disk mass of M = 2.67 • 10" 2 M Q , and a G type star (§ [3]). We 
used the 8 th ALMA configuration and a wavelength of 430 fim. A gap 
will be detected if the difference between the brightness profile and the 
auxiliary line at the local minimum of the brightness profile is three 
times larger t han the standard deviation of the noise of the observation 
(noise: see § 12.31 ). 

4. Analysis and data reduction 

We now briefly diskuss the computational analysis of the output 
from the CAS A simulator. There are three steps for the assess- 
ment of every simulated observation. First we derive the stan- 
dard deviation of a simulated observation from the correspond- 
ing noise maps (§ 12.31) . If the signal-to-noise ratio is greater than 
three, this simulated observation will be taken into account for 
further analysis. 

In the next step, the object is considered as resolved if the 
area of the synthesized beam is smaller than the object in the 
simulated observation image in a 3<x environment. This analysis 
is sufficient for the unperturbed disk model, but the other disk 
models have to fulfill a criterion that indicates the presence of 
large-scale disk structures. 

The most prominent structure a gia nt planet induces in the 
disk is a gap (iPapaloizou & Lin 1984), Because of the radial 
symmetry of the gap, it is unequivocally detectable in the radial 
brightness profile. In our setup, the disk's density models show a 
significant amount of material in the inner part of the disk (within 
the orbit of the planet). The radial brightness profile indicates 
this by a global maximum of the brightness. In the low-density 
area of the gap, the brightness profile decreases towards a local 
minimum, followed by another local maximum. We fit a line, 
parallel to the abscissa, onto this local maximum (see Fig. [6] 
for illustration). Now, a gap will be assessed as detected, if the 
difference between this auxiliary line and the local minimum is 
larger than three times the standard deviation derived from the 
noise maps (3<x environment). This order is chosen deliberately, 
because a simulated observation that does not meet a previous 
condition also does not meet all the following ones. 

The comprehensive illustration of the statistical outcome of 
the study is given in overview charts (see Fig. [J]©. They sum- 
marize the impact of the stellar type, the total disk mass, the disk 
size, and ALMA configuration on the observability in just one 
figure. If a simulated observation fulfills the criteria mentioned 
above, the mean flux per pixel of the simulated image will be 



drawn gray color-coded (in the online- version color-coded) into 
the figure. A colorbar indicates the value of the flux per pixel on 
a logarithmic scale. The presentation is only made for disks with 
masses that allow for an observation in at least one setup. 

5. Results 

We now present the main results from our studies of the ob- 
servability of planet-induced structures in circumstellar disks. 
First, the feasibility of spatially resolving the various disks and 
cut-outs from disks is diskussed. Thereafter we explore the de- 
tectability of planet-induced gaps under the influence of the in- 
dividual parameters on our setup. Based on this, we predict the 
optimal parameter combination to observe a selected object with 
ALMA. Additionally, we propose a survey of nearby, circum- 
stellar disks seen face-on. For a general overview, Figs. |7]to[8] 
show a selection of five overview charts (see § H}. In compar- 
ison, the results for an unperturbed disk are shown in Fig. |7]at 
the top. Thereafter, the outcome for a disk with a planet with 
7] = 0.001 is drawn. Until § !5.4l onlv the hydrodynamical case is 
explored. 

5.1. Feasibility of spatially resolving circumstellar disks 

Our investigations based on the unperturbed disk model (S 1 , see 
§ 12.11) outline that circumstellar disks at a distance of 140 pc can 
be spatial resolved with ALMA around all stars, disk sizes, and 
at all wavelengths in our setup. Even observations of disks with 
masses down to 2.67 • 10" 6 M o are feasible (see Fig. 13). Gen- 
erally, the number of possible observations that can resolve the 
disk with a signal-to-noise ratio greater than three increases with 
the total disk mass. In a very first approximation, the re-emission 
luminosity of a disk seen face-on depends linearly on the disk 
mass. In appendix[C]we summarize a more accurate, analytical 
prediction of the impact of the disk mass on the observability. 

Besides the number of dust particles in the disk, which de- 
pends linearly on the total disk mass, the re-emission proper- 
ties of the dust grains influence the disk luminosity. For this 
reason, we si mulat e the radiative transfer for two different dust 
species (see § 12.2b . Figure [2] shows the different absorption and 
re-emission properties of both dust populations. At a wavelength 
of 1 mm, the re-emission ability of the small dust grains is just 
1% of the ability of the large dust grains. Differences are also 
reflected in the SEDs (see Fig. [9]). While the maximum flux 
for both dust species is nearly the same, the wavelength of the 
maximum flux moves to shorter wavelengths for the smaller dust 
grains. Therefore a disk containing large dust particles re -emits 
more effectiv in the (sub)mm range, while a disk made up of 
small dust grains is more luminous in the infrared wavelength 
regime.This is also reflected by the spectral index: 



In the wavelength range between 730/mi and 3.3 mm the spec- 
tral index for the large grains is a S pec, large = 2.84 ± 0.01 and for 
the small dust grains a S pec, small = 3.96 ± 0.01. However, the flux 
of the disk in the case of smaller grain just drops to 10% of the 
larger particle case. This is because the disk mass is the same in 
both scenarios. Owing to the average volume of the dust grains 
we need more smaller dust grains in comparison to the larger 
ones to produce the same disk mass. In our study, the effective 
dust grain radius of the small particles is about 10% of the larger 
ones (flerT, small ~ 0.1 • <2 e ff, large)- Because of that, we need about 
thousand times more small dust grains for the same disk mass. 
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Fig. 7. Feasibility of spatially resolving an unperturbed disk (top) and of detecting a gap induced by a planet with a planet- to- star mass ratio of 
77 = 0.001 (bottom) depending on the disk mass, the host star, the disk size, and ALMA configuration. The disks are made up of large dust grains. 
The exposure time is two hours at a wavelength of 430 yt/m. The simulations for a combination of a total disk mass of M = 2.67 • 10~ 2 M© and a 
scaling factor k < 10 were not calculated. For the scaling factor k see Fig. \5\ for the resulting angular resolution depending on wavelength and 
ALMA configuration see Fig. [3] The scale of one pixel is 0.00064" x k. 



Although more dust in the disk influences the re-emitted flux 5.2. Detectability of planet-induced Gaps 
of the disk, the re-emission properties of the particles are more 

important within our parameter space The feasibility for a re- Thg individual parameters of our setup will influence the de- 
solved detection drops to about a half when assuming smaller tectabmt of plane t-induced gaps within the circumstellar disk, 
dust particles instead of large ones in every setup, independent In ^ foUowing we explore their impact for a hydrodynami- 
of the underlying hydro- or magnetohydrodynamical simulation. ^ disk modd wkh a ^ wim = 0Q1 Since ^ 



Article number, page 8 offTTl 



Ruge et aL: Tracing large-scale structures in circumstellar disks with ALMA 
A = 430 jum, small dust grains, planet with tj = 0.001 in the disk (SP1) 



Mass[M ] 4 7 

10 

2.67 -IO" 2 11 

19 

„ 22 
4 

O 7 

u io 

2.67- 10" 3 ^ 11 

E 19 

75 22 



Kstar 
n — i — i — i — i — r 



I I I I I I 



_i i i i i i_ 



G star 
n — i — i — i — i — r 



I I I I I I 



_i i i i i i_ 



F star 
n — i — i — r 



I I l " l I I 



_i i i i i i_ 



A star 
~i — i — i — i — i — r 



M I 'M I 



TTauri Herbig Ae 
n — i — i — i — i — r 



I I I I I I 



j i i i i i_ 



n — i — i — i — i — r 



-H- 



I I I' l I I 



W 4 8 12 16 20 24 4 8 12 16 20 24 4 8 12 16 20 24 4 8 12 16 20 24 4 8 12 16 20 24 4 8 12 16 20 24 

ALMA configuration 



n 



I 







-2 x 



— A CD 
H Q_ 

X 

-6 i 

C 
03 

-8 J 



-io § 



X - 1300 //m, large dust grains, planet with tj = 0.001 in the disk (SP1) 
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A = 330 //m, large dust grains, planet with tj = 0.001 in the disk 
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Fig. 8. Wavelength dependence of the detectability of a gap induced by a planet with 77 = 0.001 in a circumstellar disk depending on the disk 
mass, the host star, the disk size, and the ALMA configuration. Both grain species are considered. The assumed exposure time amounts to two 
hours. The simulations for a combination of a total disk mass of M = 2.67 • IO -2 M© and a scaling factor k < 10 were not calculated. For the 
scaling factor k see Fig. \5\ for the resulting angular resolution depending on wavelength and ALMA configuration see Fig. [3] The scale of one 
pixel is 0.00064" x k. 
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Fig. 9. Differences in the SED resulting from different dust species. 
The solid line represents the large dust grains and the dashed line the 
small dust grains. In both cases an unperturbed circumstellar disk of 
M = 2.67 • 10~ 2 M© is used. A G type star hosts this disk with an inner 
radius of R m = 38 AU and an outer radius of R out = 171 AU. In the 
wavelength range between 730yt/m and 3.3 mm, the spectral index for 
the large grains is ar sp ec, large = 2.84 + 0.01 and for the small dust grains 

»spec, small = 3.96 + 0.01. 
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Fig. 10. Comparison of the radial brightness distributions of identical 
disks around an A type (solid line) and a K type star (dashed line). 



tity T] respresents the planet- to- star mass ratio (see § 12. IK the 
T Tauri star (Af = 0.5 M©) is orbited by a planet of a mass of 
M = 0.5 M Jup and the Herbig Ae star (M = 2.5 M Q ) is orbited 
by a planet of a mass of M = 2.5 Mj up . 

Host star: Without other heating sources, the stellar luminosity 
and temperature define the thermal structure and the re-emission 
of our disk models. Consequently, a luminous, high tempera- 
ture star, like an A type star, provides much better conditions for 
an optimal observation of a disk in re-emission radiation than 
a K type star (for stellar parameters see § |3}. The A type star 
engages a lot more energy in the disk and heats it up more ef- 
ficiently. Because of the resulting higher dust temperature, the 
dust emission is enhanced. Therefore, the contrast between re- 
gions with and without dust increases. Figure [10] shows an ex- 
ample for the resulting radial brightness distributions of a disk 
with a planet with r] = 0.001, which in the first case is illumi- 
nated by the A type star and in the second by the K type star. 
Through the increased disk remission flux in the first case, the 



gap is detectable as a minimum in the radial brightness distri- 
bution at « 125 AU, while in the second case one only finds a 
saddle point at the location of the gap. In the parameter space 
considered in our study the gap detection rate is higher by 60% 
in the case of disks around A type stars, than in the case of a K 
type star. 

ALMA configuration: We consider a distance of 140 pc for ev- 
ery simulated disk. Therefore, the required spatial resolution of 
an ALMA configuration only depends on the size of the struc- 
ture that is intended to be observed. However, the sensitivity of 
ALMA decreases with increasing resolution. Accordingly, the 
optimal ALMA configuration is a compromise between the res- 
olution and the sensitivity. Generally the best choice for observ- 
ing a specific structure is to use an ALMA configuration with 
a synthesized beam that is slightly smaller than the structure it- 
self. Thus, it is fully resolved, and a maximum flux is caught by 
the antennas. In general, for every combination of structure size 
and its re-emitted flux, there is an optimal ALMA configuration 
for the planned observation (see Fig. [7]). For example, the vast 
majority of observations at a wavelength of 430 pm in the Fig. 
[8] is feasible with the 14 th ALMA configuration (max. baseline 
« 1600 m). This configuration is a valid compromise for a wide 
range of disk sizes. 

In this context, Fig. [TT] answers whether a structure is large 
and luminous enough to be observed by ALMA, or not. Within 
the considered model space and combining the observing wave- 
length and ALMA array configuration to the resolution, a struc- 
ture will be observable with ALMA if its combination of re- 
quired luminosity and required angular resolution is located in 
the dark area in the middle of Fig. [TT] The lower limit of this 
dark area is the detection limit of ALMA. Above, in the gray 
area, observations are possible, but these and the ones repre- 
sented through the regions on the left and right of the dark area 
are not covered by our parameter space. 

In our setup ALMA allows for detections of circumstellar 
disks down to a scale of « 2 AU (Fig. [TT] for a distance of 
140 pc). A finer resolution is limited by the sensitivity and 
position-dependent maximum exposure time. (Below an eleva- 
tion of20^_oyer the horizon an observation should not be carried 
out, IOtt&Kernl 20TTr) 



Observing wavelength and exposure time: The impact of 
wavelength and exposure time is illustrated in Fig. [12] This his- 
togram shows that it is most likely to observe a gap in a disk 
in the wavelength range from 430 pm to 2000 pm for each simu- 
lated exposure time (total wavelength range in our setup: 330 pm 
to 3300 yum). While the resolution of ALMA is higher in the 
lower part of this regime, the density structure of the disk is bet- 
ter traced at wavelengths longer than 950 pm, because of the low 
optical depth of the disk (see appendix 0), the low temperature 
(« 20 K) of the most parts of the disk, and the resulting high re- 
emission ability in the mm- wavelength regime. In the context of 
our model setup, the wavelengths 430 pm and 1 300 pm under the 
respective atmospheric conditions are good choices for search- 
ing gaps in circumstellar disks, because a maximum number of 
observations is possible (Fig. [12]). At both wavelengths, ALMA 
detects gaps down to a disk mass of 2.67 • 10" 4 M (see Fig. [7] 
bottom and [8] middle). 

In comparison to the wavelength range from 430 pm to 
2000 pm, the possibility of detecting a gap within a circumstellar 
disk decreases strongly at a wavelength of 330yum. Here, gaps 
can be traced only in the most massive disks around the most 
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Fig. 11. Overview of the correlation between resolution and sensitiv- 
ity of ALMA. The resolution combines the observing wavelength and 
the ALMA array configuration. The dark area indicates 3cr-detections 
within an exposure time of two hours. In the gray area observations 
are possible, but were not calculated in the current study. The ALMA 
detection limit is determined by the lower edge of the dark area. 



Feasible observations due to wavelength and observation time 
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Fig. 12. Overview of the influence of the exposure time h: black 
bar, 2 h: transparent bar, 8 h: gray bar) and wavelength on the feasibility 
of detecting a gap in the disk. The underlying model contains all sim- 
ulations (all disk masses, sizes, and stellar types) for a hydrodynamical 
simulation (SP1) with a planet (rj = 0.001) and large dust grains. The 
observing probability is normalized to the maximum at 430 jim and an 
exposure time of 2 h. 



luminous stars (see Fig. [8] bottom). Even with optimum ob- 
servation conditions the thermal noise dissipates the incoming 
signal too much. This effect is worse the nearer the object is lo- 
cated to the horizon. At wavelengths longer than 2000 /mi, the 
re-emitted flux of the disk decreases below the detection limit of 
ALMA (see § 15.2b or the resolution of the interferometer is no 
longer sufficient. 

The exposure time (see Fig. [12]) affects the signal-to-noise 
ratio and the coverage of the u-v plane in Fourier space. With 
the exception of the shorter wavelengths in our setup, the signal- 
to-noise ratio grows, as expected, with the square root of the 
exposure time. The decrease in the likelihood of tracing a gap 
at the shorter wavelength between a two-hour and an eight-hour 
observation is again a result of the atmosphere. During an eight- 
hour observation our object is located just « 20° over the horizon 




R [AU] 

Fig. 13. Example of a simulated observational result. Around a 
Herbig Ae star even the circumplanetary region can be detected. The 
simulated observation of the disk with M = 2.67 • 10~ 2 M© is made 
at wavelength of 430 /um with the 18 th ALMA configuration. The ex- 
posure time amounts to two hours. Dust phase: large dust grains. The 
inner radius of the disk is i? in = 14 AU, the semi-major axis of the planet 
R p i = 35 AU, and an outer radius of R out = 63 AU. 

for about one hour. Therefore, the signal-to-noise ratio decreases 
rapidly. 

In about 66% of the simulated observations, in which a gap 
was detected in a two-hour observation, this gap can already be 
traced within half an hour (see Fig. \Y2\ . An even shorter expo- 
sure time carries the risk of too low u-v coverage for the more 
extended configurations. On the other hand, an eight-hour ob- 
servation will not produce much better results than a two-hour 
one. 

Planet-disk interaction: Besides the gap, the planet-disk in- 
teraction causes structures like spiral waves and circumplane- 
tary density accumulation (Papaloiz ou & Lin| 19841) The spiral 
waves created by the Lindblad resonances dLindblad 19581) pass 
through the entire disk. They are clearly visible in the face -on 
density profiles shown in Fig. [T] However, the (sub)mm bright- 
ness contrast between the spiral arms and the background disk 
surface is too low to be indicated by ALMA at a distance of 
140 pc. But spiral arms resulting from self-gravitation of mas- 
sive, compact disks can potentially be observed by ALMA, be- 
cause of the much higher densities in this kind of disk (Cossins 
et al. 120101) . 

The circumplanetary disk is not as large in size as the gap 
or the spiral arms. However, this high densi ty area in the low 
density gap produces a high level of contrast (Wol f & D' Angelol 
2005) Therefore, it can be observed by ALMA even with ne- 
glected planetary heating (Fig. [13]). The planet as an additional 
heating source, in fact, has to be considered in a realistic model 
and thus the contrast level will even increase. The strength of 
the planet-induced d isk s tructures depends on the planet- to- star 
mass ratio r] (see § 12.11) . Besides r] = 0.001, we investigate 
the detectability of gaps for hydrodynamical simulations with 
7] = 0.0001. However, in this case planet-disk interaction can 
not be traced by ALMA. 

Generally, sufficiently massive planets perturb a circumstel- 
lar disk massively. However, a planet cannot be detected un- 
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Fig. 14. Selected examples for the detection of gap (illustrated as a bright ring at a radius of k x 5 AU) in circumstellar disks, which corresponds 
to the Butterfly Star disk model. The shown scale of the disk size has to be multiplied by the scaling factor k, shown for every line on the left. 
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ambiguously in th e SED. This has already been investigated by 
I Wolf & D'Ange lo (2005), Our results confirm this statement. 

Optimal parameters for an observation: Based on our investi- 
gations, we propose to use an exposure time of about two hours, 
a wavelength of 430 yum or 1300yum, and the 14 th ALMA con- 
figuration (max. baseline « 1600 m) for the search of planet- 
induced structures in circumstellar disks in the Taurus-Auriga 
star formation region. 

5.3. Application for real objects 

In the following we briefly diskuss the possibility of tracing 
planet-induced structures in a disk that is comparable to the 
specific, well- studied case of the Butterfly Star. The Butterfly 
Star is a perfectlv edg e-on seen Class I source in Taurus ( Wolf 
et al. 20031) For this purpose we use the mass of cut-outs from 
the Butterfly Star disk model (see appendix |Bj and disk density 
structure from a hydrodynamical disk model with a planet with 
n = 0.001. In this case the disk size already defines the disk 
mass. The masses that we consider are shown in Table [4] The 
disk size still is a function of the scaling factor k (see § [3]). Our 
results are shown in Fig. [14] With the exception of the smallest 
disks in our setup (R out = 9 AU corresponding to the scaling fac- 
tor k = 1), gaps are traced for every disk size (k = 1 is not plotted 
in Fig. [14]). The synthesized beam is about the same size in ev- 
ery image and always comparable to the corresponding gap size. 
Thus, ALMA can trace gaps with nearly all semi-major axis at 
several wavelengths with identical resolution. Besides the actual 
gap detection therefore statements on the dust in the disk are an 
additional outcome of an observation. 

We find that ALMA will be able to indicate a planet-induced 
gap in a disk at a distance of 140 pc if the planet has a semi- 
major axi s of 10AU (diameter of the gap » 2AU). Wolf & 
D'Angelo (2005) detected this structure, which was caused by 
a giant planet with a semi-major axis of 5 AU at a wavelength 
of 330 /mi at a distance of 100 pc. The main differences are the 
wavelength of the observation and the boundary conditions of 
the hydrodynamical simulations. However, for observations at a 
wavelength of 330 /mi, meaningful results are uncertain because 
of the absorption properties of the atmosphere (§ 15. 2) . In gen- 
eral, it is possible to trace planet-induced structures in 140 pc dis- 
tant circumstellar disks at AU- scale with a wavelength of about 
430 /mi and a small amount of material between star and planet, 
because of the resulting higher contrast between the gap and its 
outer edge. 

5.4. Effects of magnetohydrodynamical turbulence 

Our setup includes equivalent magnetohydrodynamical simula- 
tions for every pure hydrodynami cal simulatio n. Their classifi- 
cation and diskussion is given in lUribe et al.1 (20111) in detail. 
First we focus on a comparison between the magnetohydrody- 
namical simulation and the pure hydrod ynamical s i mulat ion of 
a disk with a planet with 77 = 0.001. lUribe et al. (2011) have 
already investigated that the presence of magnetohydrodynami- 
cal t urbulence in the disk enlarges the gap size (see also Winters 
et al. 2003 ; Nel son & Papaloiz ou 2003 ). ALMA again allows us 
to measure the size of the gap if the synthesized beam of the con- 
figuration is about half of the gap size. By this, a different gap 
size is not the result of an under-resolved observation. For ex- 
ample, Fig. [15] compares the radial brightness distribution (nor- 
malized to its maximum) of the pure hydrodynamical simulation 



Table 4. Disk dimensions and masses for the real object application. 
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Notes. Dimension and masses of the cut-outs from the Butterfly Star 
disk model (M tot = 0.07 M Q ). Besides the scaling factor k, the inner and 
outer disk radii (R in and R out ) and the planet major semi-axis (R p i an et) 
are listed. 
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Fig. 15. Measuring the gap size in the radial brightness distribution 
resulting from an ALMA observation and in the surface density distri- 
bution. The solid line represents the radial brightness profile and the 
dotted line the surface density of a disk model hosting a planet with 
77 = 0.001 without turbulence (SP1). The dashed line represents the 
radial brightness profil and the crosses the surface density of correspod- 
ing disk model with magnetohydrodynamical turbulence (SPT6). The 
central star is a Herbig Ae. 

and the magnetohydrodynamical simulation and the correspond- 
ing surface density structures (also normalized to their maxima) 
in an identical disk setup. The observation measures the gap size 
in good agreement with the surface density profiles. The differ- 
ences are mainly caused by the scale of the synthesized beam 
and the signal-to-noise ratio of the observation. 

The magnetohydrodynamical turbulence also has an impact 
in the case of simulations of planets with n = 0.0001. While 
gaps cannot be traced in the pure hydrodynamical case, it will 
be feasible to detect gaps if magnetic fields are considered. It is 
possible to detect this gap in about 10% of the cases in which a 
gap induced by a planet with n = 0.001 has already been traced. 

However, a larger gap is not a proof of the importance of tur- 
bulence and dissipative effects in the disk. It also can be formed 
by a planet containing more mass (see § [T]). Therefore, further 
observations that can determine the planet mass are needed to 
analyze the origin of the gap size, i.e. a measurement of the or- 
bital period and the semi-major axis of the planet. In general, the 
presence of magnetohydrodynamical turbulence improves the 
detectability of gaps. The observability of other planet-induced, 
large-scale structures is equivalent to the already diskussed pure 
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Fig. 16. Observability of brightness structures mimicking gap like 
features resulting from zonal flows. The disks is made of large dust 
grains around a Herbig Ae star at a wavelength of 1300//m with the 
20 th ALMA configuration. (M disk = 2.61 • 1(T 2 M Q ) 
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Fig. 17. Radial brightness profile (solid line) corresponding to Fig. 
|T6l The auxiliary line is given dotted. 

hydrodynamical case. Thus, the results of § !5.2l also apply to the 
magnetohydrodynamical case. 

Additionally, we explore the observability of zonal flows (see 
§ 12.1b in circumstellar disks. In very few cases (^20 images) 
a gap-like structure resulting from zonal flows is detectable at 
a wavelength of 1300/mi with the ALMA configurations 18 to 
22. In particular a massive and extended disk out of large dust 
particles around a bright stars are needed for a detection. Figure 
n~6l illustrates an example of these results. The corresponding 
radial brightness profile is shown in Fig. [T7] Zonal flows have 
an impact on the dust distribution in the disk and can trigger the 
formation of planetesimals (see § 12.11) . To distinguish between 
gaps that are induced by planets and gap-like features resulting 
from zonal flows, further multi wavelength observations of the 
object are necessary. 

6. Conclusion 

We have performed a large-scale parameter study of the observ- 
ability of planet-induced structures in circumstellar disks in the 



(sub)mm wavelength regime with ALMA. The basic assump- 
tions of the investigated parameter space are typical T Tauri 
objects at a distance of 140 pc and a declination of « 20°. 
Based on this, we explored the impact on the observation prop- 
erties of six central stars, nine disk sizes with outer radii from 
9AU to 225 AU, 15 total disk masses in the range between 
2.67 • 10" 7 M o and 4.10 • 10" 2 M o , different observation wave- 
lengths ([330yum, 3 mm]) and exposure times (^ h, 2h, 8h), 14 
ALMA configurations, two dust grain size distributions, and disk 
turbulence due to magnetic fields. For nearly all conceivable 
parameter combinations large-scale disk structures are detected 
by ALMA. Even without considering the planet as an additional 
heating source, the circumplanetary high density area is traced. 
Gaps are clearly identified for gap width that are larger than 
2AU. A more detail resolution of the disk structures is lim- 
ited by the sensitivity and the object's position-dependent max- 
imum exposure time. Nevertheless, ALMA allows measuring 
gap sizes and therefore distinguishing between planet- to- star ra- 
tios or underlying disk models, which both influence the gap 
size. Through the observing wavelength range of ALMA, in 
particular, large dust particles enhance the feasibility of tracing 
planet-induced disk structures. For typical protoplanetary disks, 
the disk mass affects the observability only slightly. In general, a 
wavelength of 430 jim or 1300 yum, in combination with the 14 th 
ALMA configuration (max. baseline « 1600 m) and an expo- 
sure time of about two hours, is ideal for future surveys aimed 
at detecting disk structures, in particular those induced by plan- 
ets. In magnetohydrodynamical disk models, it is more likely to 
trace planet-induced gaps, because the pre sence of a magnetic 
field enlarges the gap size (lUribe et al.l 201 11) We find that zonal 
flows are detectable with ALMA. A major fraction of the pre- 
sented simulations and more overview charts are available online 
at http : //wwwl . astrophysik . uni-kiel . de/~placid. 

To enhance the understanding of the gap and planet forma- 
tion processes and to further improve the knowledge of extra- 
solar planetary systems, a deep link of the radiative transfer to 
the hydro- and magnetohydrodynamical simulations is neces- 
sary. The radiation of the star affects the gap form and enlarges 
the vertical e xtension e specially of the outer ed ge (Jang-Condell 
& Turner 2012) , which finally influences the detectability. An- 
other goal is to implement the accurately modeled emission of 
the planet into the simulation setup. 
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Appendix A: Atacama Large (Sub)Millimeter Array 

We used the CAS A 3.2 simulator ( Petrv for the CASA Devel- 
opment Team 2012) to predict the observability of the selected 
star/planet/disk configurations. This simulation tool comes with 
28 ALMA array configuration files containig the position of 
every antenna in this configuration. In Fig. IA.1I the relation- 
ship between the maximum array baseline and the correspond- 
ing ALMA configuration is plotted. However, not only the max- 
imum baseline, but also the u-v coverage determines the quality 
of the reconstructed maps (i.e. iBuscherl 20121) For further in- 
formation about ALMA configurations and CASA see http:// 
almascience . eso . org/[ 
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Fig. A.l. Relationship between ALMA configuration and maximum 
baseline of the array. Maximum baselines derived from the CASA 
ALMA configuration files are shown as solid line. An exponential fit 
onto this data is shown as dashed line. Only the maximum baseline is 
shown, since the antennas within this maximum baseline are arranged 
in different configurations, the resolution of ALMA can also change. 



Appendix B: Deriving the model disk mass from the 
mass of the Butterfly Star disk 

We use the results IWolf et aTl (2003) applied to the modeling 
of the disk around the Butterfly Star to define the required disk 
mass in our model. The parameters derived for the disk model of 
this source are compar able to those for other circumstellar disks 
(Aftotei = 0-07 M seelMadlener et al.l 20121: iGraT e et al. 20111; 
ISauter et al. 2009 1 Scheserer et al. 2008 Wol f et al.l 20081; Wolf 
et al. 20031). IWolf et al.l (20031) used a disk density profile of the 
form 



( ^ ( R *\° I 1 
PdiskO,£) = Pol — I ex P] 2 



h(r) 



(B.l) 



where 7?* is the stellar radius and h (r) is the scale height of the 
disk: 



h(r) = hl — 



(B.2) 



From the volume integration of this disk density profile adapted 
to the Butterfly Star disk (a = 2.367, p = 1.29), we derive the 
total disk mass M(R) within the radius R: 



M(R) = 3.70 • 10" 4 M • (R [AU]) ( 



0.923 



(B.3) 



If we interpret one of our simulated and scaled disks as a cu t-out 
from the Butterfly Star disk, we can use the correlation Eq. (IB. 31) 
to estimate the total disk mass of such a cut-out. Since the outer 
radius in our studies is a function of the scaling factor k (see § 
[3]), there are nine resulting disk masses, already listed in Tableffl 

Appendix C: Optical depth and total disk mass 

The total disk mass is one of the main parameters for the observ- 
ability. Especially the highest disk masses feature good condi- 
tions for an observation (Figs. [7] and [8]). The luminosity of an 
optically thin debris disks linearlv depends on the du st mass ( Er- 
tel et al. 2012), but for optically thick disks, like protoplanetary 
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Fig. C.l. Comparison between the analytical (dashed line) and nu- 
merical model (solid line) for the optical depth in vertical direction in a 
circumstellar disk corresponding to the model for the Butterfly Star as 
diskussed in appendix ICl The disk consists of large dust grains. 
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Fig. C.2. Impact of the total disk mass on the optical depth in vertical 
direction of an unperturbed disk with an inner radius of i? in = 38 AU and 
an outer radius of R out = 171 AU at a wavelength of 1.3 mm. The disk 
consists of large dust grains. Different line styles represent different 
total disk masses. Please see the legend in the upper right of the figure. 



disks, this estimation is more difficult. For the optical thickness 
r y along a way dz in the vertical direction of the disk with a 
mass density of the dust pd us t and a disk density profile Pdisk, we 
achieve 



J n ■ Gext(v) • a\ s 

— oo 

+oo 

3 Gext(v) 



Pdisk 



4 e ff Pdust 



Pdisk dz = 



Pdust ' (f^eff) 
3 Gext(v) 



dz 



4 e ff Pdust 



(C.l) 



(C.2) 



Here I is the surface density of the disk and <2 e xt the extinction 
efficiency. Based on this analytical estimation we test a numer- 
ical implementation in MC3D for the disk around the Butterfly 
Star (see Fig. IC.ll) . For the density structures from the hydro- 
and magnetohydrodynamical simulations (see § 12. IK we find 
that the optical thickness of the disk in the vertical direction is 
too low to influence the radiative transfer in the Rayleigh-Jeans 
wavelength regime (see Fig. IC.2I) . 
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Table C.l. Values and standard deviation for the fitted y-parameter from 
Eq.lC4l 



total disk mass 



2.67 


10" 


■ 2 M 


0.64 ± 3.96 • 


10" 


-5 


2.67 


10" 


• 3 M 


0.58 + 3.56- 


10" 


-5 


2.67 


10" 


" 4 M 


0.49+ 1.02- 


10" 


-5 


2.67 


10" 


■ 5 M 


0.43 + 1.53- 


10" 


-7 


2.67 


10" 


■ 6 M 


0.41 + 8.63 • 


10" 


-8 


2.67 


10" 


■ 7 M 


0.42+ 1.33- 


10" 


-7 



L0° 



Rayleigh-Jeans-approximation 









































































Model with fitting errors 
Linear relation 
Relation from SEDs 











1Q -7 1Q -b 1Q -5 1Q -4 lQ -6 1Q -2 1Q -1 

Mass [M ] 

Fig. C.3. Influence of the total disk mass on the re-emitted flux of an 
unperturbed circumstellar disk. A comparison between a linear model 
(dotted line) that is useful for very low mass disks, a model on the ba- 
sis of the Rayleigh- Jeans-law (solid line), and the flux simulated in the 
radiative transfer (dashed line) is shown. Every disk has an inner radius 
of Ri n = 38 AU and an outer radius of R out = 171 AU, and the hosting 
star is a G type star (§ [3]). The disk consists of large dust grains. 



With a constant disk size, dust grain size, and dust mass den- 
sity, the only process for increasing the mass within the disk is 
to increase the number of particles (linearly dependence). While 
the disk is optically thin in the vertical direction (see Fig. IC.2b , 
the effective luminous surface of the disk A depends on the num- 
ber of dust particles within the model space, so we assume 



A=A 



o • 



M 
Mo' 



(C.3) 



This fact is sufficient for debris disks, but for a circumstellar 
disk the optical thickness along a line of sight from the star has to 
be considered, too. Besides the linear increasing effective lumi- 
nous surface of the disk, we find that the temperature profile of 
the disk midplane can be used to approximate the re-emitted flux 
in the Rayleigh- Jeans-regime. We assume the following temper- 
ature profile: 



T(R) = C R- y . 



(C.4) 



We take Co = 300 K and ye [0, 1] (y depends on the dust mass 
M). Using the Rayleigh- Jeans-law: 



BAT) 



I 4 



(C.5) 



and the gray-body approximation of the dust, we get the flux F 
of a radial symmetric disk seen face-on by inserting Eq. (IC.4I) 



into Eq. dC.5b and integrating over the radius coordinate: 



2 7T 2 C &b 

A 4 



Rom 

GabsW'CoA J R~ y dR. 



(C.6) 



Neglecting the case y = 1, we can estimate the flux F: 



F = 



2 7T 2 C &b 

A 4 M 



GabsW)-A -Co 



l-y 
out 



R 



l-y\ 



1-7 



• M. 



(C.7) 



This result is comparable to lBeckwith et al.l (20001) who used a 
constant disk temperature. For testing our model, we use six cir- 
cumstellar disks of the same size but different disk masses. From 
the SEDs we get their flux at a wavelength of 1300yum. There- 
after we fit the parameter y from Eq. (IC.4l) to the corresponding 
temperature profiles. The results are shown in Fig. IC.31 and the 
fitting values of y are listed to the left of the figure. The benefit of 
the method is the easily prediction of the disk temperature profile 
and its mass from a spatially resolved millimeter observation. 
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